GLOBAL WELL-POSEDNESS FOR THE EULER-BOUSSINESQ SYSTEM 

WITH AXISYMMETRIC DATA 



TAOUFIK HMIDI AND FREDERIC ROUSSET 



Abstract. In this paper we prove the global well-posedness for the three-dimensional Euler- 
Boussinesq system with axisymmetric initial data without swirl. This system couples the Euler 
equation with a transport-diffusion equation governing the temperature. 
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1. Introduction 

Boussinesq systems are widely used to model the dynamics of the ocean or the atmosphere. They 
arise from the density dependent fluid equations by using the so-called Boussinesq approximation 
which consists in neglecting the density dependence in all the terms but the one involving the 
gravity. This approximation can be justified from compressible fluid equations by a simultaneous 
low Mach number/Froude number limit, we refer to [H] for a rigorous justification. In this paper 
we shall assume that the fluid is inviscid but heat-conducting and hence the system reads 



(1) 



dtv + V ■ Vv + Vp = pcz, (t, x) G 
dtp + v -Vp- Ap = 
divv = 0, 

V\t=0 = Vo, p\t=o = po. 
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Here, the velocity v = {v^,v'^,v^) is a three-component vector field with zero divergence, the scalar 
function p denotes the density or the temperature and p the pressure of the fluid. Note that we 
have assumed that the heat conductivity coefficient is one, one can always reduce the problem to 
this situation by a change of scale (as soon as the fluid is assumed to be heat conducting) which 
is not important for global well-posedness issues with data of arbitrary size that we shall consider 
here. The term pcz where = (0,0,1)* takes into account the influence of the gravity and the 
stratification on the motion of the fluid. Note that when the initial density po is identically zero 
(or constant) then the above system reduces to the classical incompressible Euler equation: 

{dtv + V ■ Vv + Vp = 
divt; = 
V\t=0 = Vq. 

From this observation, one cannot expect to have a better theory for the Boussinesq system than 
for the Euler equation. For the Euler equation, a well-known criterion for the existence of global 
smooth solution is the Beale-Kato-Majda criterion [3]. It states that the control of the vorticity of 
the fluid w = curlu in Lj^^^(M+, L°°) is sufficient to get global well-posedness. In space dimension 
two, the vorticity uj can be identified to a scalar function which solves the transport equation 

dfUJ + V ■ Vuj = 0. 

From this transport equation, one immediately gets that 

||^^(i)||Lp < II^^oIIlp 

for every p > 1 and hence the global well-posedness follows from the Beale-Kato-Majda criterion. 
In a similar way, the global well-posedness for two-dimensional Boussinesq systems which has 
recently drawn a lot of attention seems to be in a satisfactory state. More precisely global well- 
posedness has been shown in various function spaces and for different viscosities, we refer for 
example toPElElIIIlIiailelllTlIIHlIiniEI]. In particular, for the model ([I]) in 2D, the main idea 
is that by studying carefully the coupling between the two equations and by using the smoothing 
effect of the second equation, it is still possible to get an a priori estimate in L°° (or in i?^ ^ ) for 
u) and hence the global well-posedness. 

In the three-dimensional case, very few is known: even for the Euler equation, the vorticity uj solves 
the equation 

(3) dttd + V ■ Vu = UJ ■ Vv 

and the way to control the vortex stretching term uj ■ Vv in the right-hand side is a widely open 
problem. Nevertheless, a classical situation where one can get global existence is the case that v is 
axisymmetric without swirl [29j, [22j . Our aim here is to study how this classical global existence 
result for axisymmetric data for the Euler equation can be extended to the Boussinesq system ([1]). 
Before stating our main result, let us recall the main ingredient in the global existence proof for the 
Euler equation with axisymmetric data. The assumption that the vector field v is axisymmetric 
without swirl means that it has the form: 

(4) v{t,x) = v^{t,r,z)er +v^{t,r,z)ez, x = {xi,X2,z), r = (xf + xl)^ , 

where (^er,eg,ez) is the local basis of corresponding to cylindrical coordinates. Note that we 
assume that the velocity is invariant by rotation around the vertical axis (axisymmetric flow) and 
that the angular component v^ of v is identically zero (without swirl). For these flows, the vorticity 
is under the form 

UJ = {d^v"" - drv'')eg := ujgee 
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and the vortex stretching term reads 



OJ • \/v = UJ. 

r 



In particular ujg satisfies the equation 

(5) dtu;e + v ■ Vujq = —^e- 

r 

The crucial fact is then that the quantity C •= solves the transport equation 
from which we get that for every p € [1, oo] 

(6) \m\\Lv < iicoIIlp. 

It was shown by M. Ukhoviskii and V. Yudovich [29j and independently by O. A. Ladyzhenskaya 



p2j that these new a priori estimates are strong enough to prevent the formation of singularities in 
finite time for axisymmetric flows without swirl. More precisely global existence and uniqueness was 
established for axisymmetric initial data with finite energy and satisfying in addition ujq ^ L'^n L°° 
and ^ £ n -L°°. In terms of Sobolev regularity these assumptions are satisfied if the velocity vq 
belongs to -ff* with s > |. This condition was improved more recently. In |28j, it was proven that 
global well-posedness still holds if vq is in H'^ with s > | (note that this is the natural regularity 
requirement for the initial velocity in the Sobolev scale in view of the standard local existence 
result) and in the recent work [10], Danchin has obtained global existence and uniqueness for initial 
data such that G L^'^ n L°° and Co £ L^'^(here, L^'^ denotes the Lorentz space, the definition 
of the Lorentz spaces U'''' as interpolation spaces is recalled below). Their proof is based on the 
observation that one can deduce from the Biot-Savart law, the pointwise estimate 

v"" 1 

By convolution laws in Lorentz spaces (again recalled below) and ([U]), this yields the estimate 

-it) ^ <IIC(t)llL3.i<||Co||L3,i. 

Since one gets from a crude estimate on ([5]) that 

\\Mt)\\L^ < ||a.o(t)||L-e/olK/HlLoo^ 

the global well-posedness in s > 5/2 (the assumption (^q G L^'^ is automatically satisfied) then 
follows from the Beale-Kato-Majda criterion. It is actually possible, as shown in [2], to get global 

well-posedness in the critical Besov regularity, that is, vq G B^^ , \/p € [l,oo], in the sense that it 
is possible to propagate globally the critical Besov regularity if (q € L^'^. 

Our aim here is to extend these global well-posedness results to the Boussinesq system ([T|). Our 
main result reads: 

Theorem 1.1. Consider the Boussinesq system ([1]). Let s > |, fo G -ff** be an axisymmetric 
divergence free vector field without swirl and let po be an axisymmetric function belonging to H'^^'^n 
with m > 6 and such that r^po G L^. Then there is a unique global solution {v,p) such that 

{v,p)GCiR+;H') X (c{R+;H^-^nL"^)nLl^{R+;W'n) 

and r'^p e C{R+;L'^). 



Let us give a few comments about our result. 
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Remark 1.2. By axisymmetric scalar function we mean again a function that depends only on the 
variables (r, z) but not on the angle 9 in cylindrical coordinates. One can easily check that for 
smooth local solutions, if (poi^o) is axisymmetric (and vq without swirl), this property is preserved 
by the evolution. 

Remark 1.3. The assumption on the moment of p is probably technical. The control of the moments 
of p are needed in our proof in some commutator estimates (see ([9]) for example). 

Note that in view of the proof for the Euler equation, the crucial part is to get an a priori estimate 
for in L^'^. The equation for = ujg/r becomes 

(8) dtC + v-VC = -^ 

r 

and consequently, the main difficulty is to find some strong a priori estimates on p to control the 
term in the right-hand side of ([8|) . The rough idea is that on the axis r = the singularity - scales 
as a derivative and hence that the forcing term drp/r can be thought as a Laplacian of p and thus 
one may try to use smoothing effects to control it. We observe that if we neglect for the moment the 
advection term v-'Vp in the equation of the density then using the maximal smoothing effects of the 
heat semigroup we can gain two derivatives by integrating in time which is exactly what we need. 
From this point of view we see that our model is in some sense critical for the global well-posedness 
analysis. The main difficulty if one wants to use this argument is to deal with the advection term. 
Indeed, the only control on v that we have at our disposal is a L'^^L'^ estimate (which comes from 
the basic energy estimate) and this is not sufficient to obtain an estimate for D^p in Lj^^{LP) by 
considering the convection term as a source term and by using the maximal smoothing effect of the 
heat equation. Even more refined maximal regularity estimates on convection-diffusion equations 
( [H] ) IE] for example) do not seem to provide useful information when the control of the velocity 
field is so poor. Consequently, our strategy for the proof will be to use more carefully the structure 
of the coupling between the two equations of ([1]) in order to find suitable a priori estimates for 
(C, p)- Since the coupling between the two equations does not make the original Boussinesq system 
well suited for a priori estimates, our main idea is to use an approach that was successfully used 
for the study of two-dimensional systems with a critical dissipation, see \18\ [T9] and the Navier- 
Stokes-Boussinesq system with axisymmetric data [20]. It consists in diagonalizing the linear part 
of the system satisfied by C and p. We introduce a new unknown T which here formally reads 

r = c + ^A-V 

r 

and we study the system satisfied by (F, p) which is given by: 

dtT + V ■ S7T = - [^A-^ ,v V]p, dtp + v-Vp = Ap 
where [^A~-^,f • V] is the commutator defined by 

[^A-\v V]p= ^A-\v ■Vp)-v V(^A-V). 

Note that if we forget the commutator for a while, we immediately get an a priori estimate for 
r for every p from which we can hope to get an estimate for if the operator ^A~^ behaves 
well. 

To make this argument rigorous, we need first to study the action of the operator ^A~^ over 
axisymmetric functions. This is done in Proposition ()2.9p where we prove that this operator takes 
the form ^ 

-yA-^ = Y,aij{x)nij 

hi 
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where TZij = didjA ^ are Riesz operators and the functions aij are bounded. This yields that 
^A~^ acts continuously on L^'^ and hence that 

||^A-V(t)|L3,i<l|p(t)llL3,i<||po||L3.i. 

It follows that the control of F is equivalent to the control of C in L^'^. Now it remains to estimate 
in a suitable way the commutator term [^A~^, u • V]p which is the main technical part. It seems 
that there is no hope to bound the commutator without using unknown quantities because there 
is no other known a priori estimates of the velocity except that given by energy estimate which is 
not strong enough. We shall prove (Theorem 13. ip that 

(9) [(.dr/r)A-'^,vV]p <\\uJe/r\\L3,i{\\pxh\\B0^nL2 + \\p\\ i )• 

This estimate is the heart of our argument, its proof combines the use of paradifferential calculus 
and some harmonic analysis results and also requires a careful use of the property that velocity v 
is axisymmetric without swirl in the Biot-Savart law. 

The main reason for which we need some moments of p in the right-hand side of ([9]) is that we 
want an estimate of the commutator invoving ujg/r and not uj. 

In the right hand side of ([9]), ||/?|| i and ||px/j||^2 can be controlled in terms of the initial data 

^2.1 

only by using the smoothing effect of the convection-diffusion equation for p and standard energy 
estimates. Consequently, from this commutator estimate, we obtain that 

(10) ||C(t)||L3,i <C(t)e 

and the next difficult step is to control ||/ox/i||j;^i^o ^. This is done in two steps. The first step is 
to get a global estimate of pxh in terms of the initial data only and then in a second step, we 
shall prove a logarithmic estimate for the ^ norm of XhP in terms of the L^'^ norm of 
For the first step, let us observe that / = pxh solves the equation 

(11) dtf + v-Vf - Af = v^p-2VhP. 

Note that for the moment, we only have at our disposal the standard energy estimate for v (thus we 
control ||f 11^00^2 only), consequently to obtain an L°° estimate for / we need to use an 
estimate for the convection-diffusion equation since the source term in the right-hand side can be 
estimated only in L^. Note that the convection term cannot be neglected (again because of the 
weak control on v that we have at this stage) and hence this estimate cannot be obtained from 
heat kernel estimates. We shall obtain this estimate by using the Nash-Moser-De Giorgi iterations 
[E]) |1S]) [21] ■ Indeed, the main interest of this approach is that since it is based on energy type 
estimates, the convection term does not contribute. A general result is recalled in the Appendix 
\K[ For technical reasons, some higher order moment estimates which are easier to obtain are also 
needed, they are stated in Proposition 14.21 

Once the estimate of ||/OX/i||i,oo is known in terms of the initial data, one can establish logarithmic 
Besov space estimates for the convection-diffusion equation ([TT]) by using a special time dependent 
frequency cut-off of x^p where we combine the estimate with some smoothing effects for x^p. 
This yields (see ([56]) ) 



(12) llpx/^ILi^o^ ^ <Co(t)(l + ^*/i(r)log(2 + ||CkooL3.i)dT 

where Co(t) is a given continuous function and h is some L\^J^j^) function. We point out that the 
use of the moment of order two is due to the treatment of the commutator [A|j,t; • V]{xhp) 

which appears when we deal with the smoothing effects. 
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The combination of the estimates (jl2p and (jlOp with Gronwall inequahty ahows to control ||C(0IIl3,i 
globally in time. 

The final step is to deduce, as for the incompressible Euler equation, from the control of ||C||l^l3.i 
an estimate of ||a;||Lj«ioo and of HVvHi^ioo. This is the aim of Proposition B3] and Estimates 
in Sobolev spaces then follow in a rather classical way. 

Once a priori estimates for sufficiently smooth functions are known, the result of Theorem 11.11 
follows from an approximation argument. 

The paper is organized as follows. In section [2] we fix the notations, give the definitions of the 
functional spaces, in particular Besov and Lorentz spaces, that we shall use and state some of their 
useful properties. We also study the operator ^A~^ in Proposition 12. 9i Next, in section [3l we 
study the commutator [^A~^,t; • V] . In section HI we turn to the proof of a priori estimates for 
sufficiently smooth solutions of ([l]). We first prove in Proposition 14. II some basic energy estimates, 
next, we study the moments of p in Proposition 14.21 and then we control HCIIl^.i in Proposition 
14.41 Lipschitz and Sobolev estimates are finally obtained in Proposition 14.51 and Proposition 14.61 
In section [5l we give the proof of Theorem 11.11 we obtain the existence part by using the a priori 
estimates and an approximation argument and then we prove the uniqueness part. Finally, the 
appendix |A] is devoted to the proof of a priori estimates for convection-diffusion equations by the 
Nash-De Giorgi iterations which are needed in the estimate of the moments of p. In the appendix 
|B]we give the proof of Lemma |2 . 71 which is a technical commutator lemma used in several places. 

2. Preliminaries 

2.1. Dyadic decomposition and functional spaces. Throughout this paper, C stands for some 
real positive constant which may be different in each occurrence and Cq denotes a positive number 
depending on the initial data only. We shall sometimes alternatively use the notation X for 
the inequality X < CY. 

When B is a Banach space, we shall use the shorthand Lj,{B) for LP(0,T, B). 

Now to introduce Besov spaces which are a generalization of Sobolev spaces we need to recall the 

dyadic decomposition of the unity in the whole space (see [7] ) . 

Proposition 2.1. There exist two positive radial functions x G ^?(M'^) and tp £ 'D(W^\{0}) such 
that 

(1) X(6 + E '^(2-^0 = 1, ^ < X\0 + E < 1 G M3, 

qeN qeN 

(2) supp (p{2'P-) n supp p{2-i-) = 0, if \p-q\> 2, 

(3) q>l^ suppx n supp 9?(2"'?) = 0. 

For every u G 5'(]R'^) we define the nonhomogeneous Littlewood-Paley operators by, 

A^iu = x(D)u; \/q G N, AgU = (p{2''^D)u and SgU = ^ ^i^- 

-i<j<9-i 

One can easily prove that for every tempered distribution u, 

(13) n = E N^- 

<?>-! 

In the sequel we will frequently use Bernstein inequalities (see for example [7]). 

Lemma 2.2. There exists a constant C such that for /cGN, l<a<6 and u G L^, we have 

sup Wd'^SguW Lb < C^'2'?(^'+3(^-i))||5^u||La, 

\a\=k 
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and for g E N 



C-^2i^\\/\qu\\L'^ < sup ||5"Agu||La < C^2'''=||Agn||L- 

\a\=k 



The basic tool of the paradifferential calculus is Bony's decomposition [4J. It distinguishes in a 
product uv three parts as follows: 

uv = TuV + TyU + Tl{u, v), 

where 

TuV = Sq-luAqV, and Tl{u, v) = AqUAqV, 

q q 

1 

with Aq = ^ Aq+i. 

i=-l 

The term TuV is called the paraproduct of by u and TZ{u, v) the remainder term. The main interest 
of the paraproduct term is that each term Sq-iuAqV has the support of its Fourier transform still 
localized in an annulus of size 2'^ and thus T^v is a sum of almost orthogonal functions. 
Let {p, r) £ [1, +oo]^ and s G M, then the nonhomogeneous Besov space Bp ^. is the set of tempered 
distributions u such that 

:= (2'"\\Aqu\\Lp) < +00. 

We remark that the Sobolev space coincides with the Besov space -622- Also, by using the 
Bernstein inequalities we get easily the embeddings 

Bp^,ri ^ Bp^,r2''' , Pi<P2 and n < r2. 

Finally, let us notice that we can also characterize spaces in terms of the dyadic decomposition, 
see [27j. For p +00 [, there exists C > such that: / belongs to LP if and only if (Aq/)q>_i E 
and 



q>-l q>-l 



LP 



where 1 < Po < P < Pi ^ 00, 6 satisfies ^ = + ^ and 1 < g < 00. 



2.2. Lorentz spaces and interpolation. For p g]1,oo[, g E [l,+oo], the Lorentz space L^'^ can 
be defined by real interpolation from Lebesgue spaces: 

(LP«,LPi)(,,,) = Lf''', 
1 _ 1- 

From this definition, we get: 

(15) LP'" ^ LP'''' , LP'P = LP 
for every 1 < p < 00, 1 < q < q' < 00. 

Lorentz spaces will arise in a natural way in our problem because of the following classical convo- 
lution results, for the proof see for instance \23\ 126] . 

Theorem 2.3. For every a, < a < d, pi e]1, +oo[, qi E [1, +00], such that 1 + ^ = ^ + ^ and 
— = — + —, there exists C > such that 

<?1 92 (?3 ' 

(16) \\f * QHlPI'^ < C\\f\\LP2>12 \\g\\LP3'13. 

Moveover, in the case that pi = <x, we have 

(17) ||/*5||l-(r<^)<C||/||^^,^^^^^||<7||^_^,, 
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In particular, by using this result and the fact that belongs to L2'°°(M^), we have that 

(18) ||VA-V||l-(ir3) < ||/||i3,i( 

and thanks to the pointwise estimate ([7]) that 



r 



, <IICIIl3,i. 



To establish some functional inequalities involving Lorentz spaces the following classical interpola- 
tion result (see [23] for example) will be very useful. 

Theorem 2.4. 

Let 1 < pi < P2 < oo, 1 < ri < r2 < oo,q £ [1, oo] and T be a linear bounded operator from L^' to 
U'K Let g]0, 1[ and p,r such that | = ^ + '^'^^ r ~ n" ^ ^^^'^ ^^•^'^ bounded from 

LP'i to L'^'i with 

As a consequence, we obtain the following results. 

Proposition 2.5. For 1 < p < +oo, q G [l,+oo], i/ien exists a constant C > such that the 
following estimates hold true 

(1) ||tiv||LP,q < C||n||Loo llvll/^p.q, 

(2) ||T„t;||iP,9 < C||n||icx; 

(3) Let us define the Riesz transform IZij = didjA~^,i,j £ {1,2}, then 

\\Ti'ijU\\LP,q < C||ti||LP,iJ . 

(4) For s > ^ we have ^ L^'^. 
Proof. 

(1) For a fixed function u £ L°°, the linear operator T : v ^ uv belongs to C{LJ',LP) with norm 
smaller that ||ti||L°° and hence the result follows by interpolation from Theorem 12.41 

(3) In a similar way, for every p g]1,+oo[, TZij £ C{LP,LP) thanks to the Calderon-Zygmund 
theorem and hence (3) follows again by using Theorem 12.41 

(2) To establish the inequality, it is again sufficient thanks to Theorem 12.41 to prove that for 
u £ L°°,v £ LP we have ||r„f||LP < C'||^*||l°° I^Ulp- For this last purpose we will make use of the 
maximal functions tool. We will start with some classical results in this subject. For a locally 
integrable function / : — )• M, we shall define its maximal function Mf hy 



Mf{x) =sup4 / \f{y)\dy. 

r>0 JB(x.r) 



B{x,r) 

From the definition we get 

(20) 0<M{fg){x)<\\g\\L^Mf{x). 

It is well-known that A4 maps continuously L^ to itself for p g]1,oo]. Moreover, we have the 
following lemma. We refer to pY] for a proof. 

Lemma 2.6. (1) Let ijj £ 5(M^) and define = e~^ip{e'^x) for e > 0. Then there exists 

C > such that for every p £ [l,oo], e £ (0, 1], we have 

sup\^|Je*fix)\<CMfix). 

In particular, we have 

sup \Agf{x)\ < CMf{x). 

q>-l 



(2) Let p g]1,oo] and {fq,q > —1} be a sequence belonging to L^l'^. Then we have 



X 



1 

|2\ 2 



< 



LP 



c||(EI/. 



1 

|2\ 2 



g q 
Let us now come back to the proof of (2). By using ()14p . we have 



LP 



1 

2\ 2 



\\TuV\\lp < \\[Y1 

i>-i 

This yields according to Lemma 12.61 and ([20 



LP 



i>-i li-<?l<4 



LP 



\TuV\\lp 

^ IK E ( E Ms,~iuA,v: 

i>-i |j-g|<4 



LP 



< 



\u\\l° 



2 



LP 



9>-l 



2\ 2 



< 



LP 



u\\l°°\\ f ^ (Aq 



1 

,2\ 2 



5>-l 



LP 



< 



rs. ,|^i||L°°|P||LP 

where the last estimate follows from a new use of (jl4p . This ends the proof of (2). 

(4) This embedding follows from Sobolev embeddings combined with Theorem 12.41 This is left to 

the reader. □ 

2.3. Some useful commutator estimates. This section is devoted to the study of some basic 
commutators which will be needed in our main commutator estimates, especially in Theorem 13.11 
and Proposition 13.21 . Our first result reads as follows. The proof is postponed to Appendix [Bl 



Lemma 2.7. Given {p,r,p,m) G [l,+oo]^ such that 

1111 



1, 



Lemma 2.8. Let p,m, p G [l,+oo] such that = ^ + j;- Then, there exists C > such that for 



1 -\ — = 1 1 — , p > r and p > 3(1 ). 

p m p r r 

Let f,g and h be three functions such that V/ G LP,g £ and xJ-^^h G L^. Then 

\\[h{B),f]g\\^^<C\\xT-'h\\Lr\\Vf\\Lp\\g\\Lr^. 

where C is a constant. 

As an application of Lemma 12.71 we get the following commutator estimates. 

1 _ 1 I 1 

V/ £ LP,g £ and for every q G NU {-1} 

\\[Ag,f]g\\^,,,<C\\Vf\\Lp\\g\\L-, 
with the following definition Hv^H^i.p = ||V(^||2,p. 
Proof. We write for i = 1,2, 3, 

di{[A„ f]g) = [d,A„ f]g - difA.g = [h,{B), f]g - OifA^g, 
with hq{C) = 2''(/>(2~«^), and G ^(M^). Using Lemma [221 we get 

\\[hq{B),f]g\\LP < C\\xT-^h,hi\\Vf\\Lp\\g\\L^<C\\Vf\\Lp\\g\\L^ 
For the other term, the Holder inequality yields 

\\difAgg\\LP < C\\Vf\\Lp\\Agg\\Ln. 

< C||V/||lp||5||l-. 



□ 



2.4. Some algebraic identities. We intend in this paragraph to describe first the action of the 
operator —A^^u over axisymmetric functions. We will show that it behaves like Riesz transforms. 
The second part is concerned with the study of some algebraic identities involving some multipliers 
which will appear in a natural way when try to study our main commutator ^dr/r)A^^,v • V]/9. 

Proposition 2.9. We have for every axisymmetric smooth scalar function u 

2 2 

(21) {dr/r)A-\{x) = ^Unuix) + ^7e22n(x) - 2^7^l2n(x), 

with TZij = dijA~^. Moreover, for p ^]l,Qo[,q G [l,oo] there exists C > such that 

(22) \\{dr/r)A-\\\LP.', < C\\u\\lp,i. 

Proof. We set / = A^^u, then we can show from Biot-Savart law that / is also axisymmetric. 
Hence we get by using polar coordinates that 

(23) duf + d22f = {dr/r)f + drrf 

where 

dr = —di + —02- 
By using this expression of dr, we obtain 

drr = {^^^ + ^82)' = 9.(^)91 + dr{^)d2 + 4^11 + ^^22 + ^^i^- 
„2 ^2 



-^0\\ + ^^22 H 5— ai2 



Since 



,x 



9^(^) = o, Vie {1,2}. 



r 

This yields by using (p3]) that 

^/ = (1 - 4)^11/ + (1- 4)^22/ -^9i2/ 

= 4^11/ + 4^22/ -^ai2/. 

To get ([2T|) . it suffices replace / by A^^u. 

The estimate (|22|) is a consequence of (|21|) and the estimates (1) and (3) of Proposition (|2.5|) since 
for every i,j e {1,2}, ^ e L°°. □ 

We shall also need the following identities and estimates. 

Lemma 2.10. For every f G 5(M^,M), we have 

(1) Fori, J G {1,2,3} 

A-\xidjf)=XidjA-^f + Cijf 
where Cijf = —2TZijA~^f. Moreover, we have the estimates: 

(24) ||VA,/||l- <C||/||l3,i, 

(25) llV^A.-ZlliP,, <C||/||iP,,, pe]i,+ooUG [l,+oo] 

(2) Fori,j,/cG {1,2,3} 

'R'ijixkf) = XkTZijf + dljf, 

with 

4 := -2dkA-^nij + 5^kdjA-^ + 5jkdiA-^ 
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where 5ij denotes the Kronecker symbol. Moreover we have the estimates 

(26) ||4/||i^ <C||/||i3.i, 

(27) I|v4/||lp.. < C\\f\\Lv,„ p e]l, +oo[, q e [l, +00]. 
Proof. (1) We first expand 

A(x,a,A-V - 27^,,•/) = 2dijA-^f + x^djf-2n^jf 

= Xjdjf 

This yields 

A-\xjdjf) = XidjA-^f - ITZijf + P{x), 
witli P a harmonic polynomial. We can easily see that the r.h.s of this identity and TZijf are 
decreasing at infinity. Thus to prove that P is zero it suffices to prove that XidjA~^f goes to zero 
at infinity. Since 

\xidjA f\ < \xj\ / I 72dy< / I 7dy+ / T^dy. 

Using Proposition 12.31 we get that xidjA^^ f G L^, for every p > 3. Hence we get P = 0. 

The estimates ()24p , (125p are a direct consequence of the above expression and (I18[) and the estimate 

(3) of Proposition 12. 5i 

(2) We use the same idea as previously. We first get the identity 

A4/ = A(7^jJ(xfc/) - XkTZijf) = dijixkf) - 2dkTlijf - Xkdijf 
= Sikdjf + 6jkdif - 2dkUijf 
and by the same argument as above, we finally obtain that 

Tlijixkf) - XkTZijf = 6ikdjA~^f + 6jkdiA~^f - 2dkTZijA~^f. 

The estimates ([26j) , (p7|) are a direct consequence of the above expression and (fT8|) and the estimate 

(3) of Proposition ESI □ 

3. Commutator estimates 

3.1. The commutator betw^een the advection operator and ^A^. hi this part we discuss 
the commutation between the operators ^A^^ and v • V. This is a crucial estimate in order to get 
better a priori estimates for the solution of ([T]) by using our transformation. Our result reads as 
follows. 

Theorem 3.1. Let v be an axisymmetric smooth and divergence free without swirl vector field and 
p an axisymmetric smooth scalar function. Then we have, with the notation Xh = (xi,X2), that 

[{dr/r)A-'^,vV]p < WuJe/rWL'i^iWpxhWBO^ ^nL2 + IIpII 1 )• 
^ ' ' ^2,1 

Proof. Since the functions p and v ■ Vp are axisymmetric then using the identity of Proposition [27 
we have 

9 2 2 

{dr/r)A^^p{x) = -|7^llp(x) + ^7^22/o(x) - 2 — —TZupix) := aij{x)TZijp{x) 



and also 



{dr/r)A'\v ■ Vp)ix) = a,j{x)'R,,j{v ■ Vp){x). 
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Since v has no swirl and the functions aij do not depend on r and z, we have for every 1 < i,i < 2 

V ■ Vaij{x) = v^drttij + v^d^ttij = 0. 
Consequently our commutator can be rewritten as 

2 2 

[{dr/r)A~'^,vV]p{x) = ^ aij{x)[nij,v -Vjp = ^ aij{x)div{[TZij,v]p} 

i,j=l i,j=l 

where we have used the fact that v is divergence free to get the last equality. By using that ajj G L°° 
and the estimate (1) of Proposition 12. 5| we first obtain that 



L3,l 



2 

(28) \\[{dr/r)A-\vV]p\\^^^^ < \\div[[n,„v]p 

The terms di[[Tlij,v^]p) and d2{[Tlij,v'^]p) can be treated in same way and hence, we shall prove 
the estimate of the first one only. The estimate of c^s ([T^jj, which is easier will be done in a 
second step. 

• Estimate o/ 5i ([T^jj, . Since v is divergence free, we have that Av = — V A w. Hence for 
axisymmetric flows (where in particular oj = uje^e), we obtain that 

v\x) = A-^dsu^ = A~^ds{xi{ue/r)). 

Applying Lemma l2.im -(l) we get 

(29) v^ix) = xiA-^dsioje/r) + Ciue/r), with £ = -Idi^A''^ 

(we omit the subscript ij for notational convenience). Consequently the commutator can be rewrit- 
ten under the form 

di{[nij,v^]p} = al([7^iJ,£(u;e/r)]p) +al([7^iJ, {A-^ds{u:e/r))xi\p^ 
= 5l([7^iJ,£(^^e/r)]p) + al([7^iJ, {A-^d3{uje/r))]xip^ + di(^{A-^d3iuje/r)) 
= di({A-^d3{u;e/r))C}jp) +di([nijX{coe/r)]p) +di([nij, {A-^ds{iOe/r))]xip) 



(30) =1 + 11 + 111 

where we have used the identity (2) of Lemma l2.10[ 
Estimate ofl. We write 



(31) di (dsA-^iJe/r) C\^p) = TZniue/r) C\^p + d:iA-\uje/r) d^L\p. 
By using (1) and (3) of Proposition 12.51 and (j26p we have 

pi\?,{ujelr) C\jp\iz,i < ||7^l3(we/r)||L3,l < C||a;e/r||L3,i ||p||i3,i 

and by using Proposition I2.5I -(1) and (fT8|) . (p7|) . we also obtain 

\\d3A-^{uj0/r)diCljp\\L3,i < Wd^A'^ {uj0 / r)\\ \\diC}jp\\L3,i < C||tJe A||l3,i ||p||l3,i . 
Combining these estimates we find 

(32) ||I||l3,i < C\\uje/r\\L3,i\\p\\L3,i. 
Estimate of II. We will use Bony decomposition 

II = IIi+Il2+Il3, 
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with 



III = 5i J];[7^^J, Vi('C(^eA))]Agp 



q>Q 



g>0 

II3 = di^[TZij,Ag{C{oJe/r))]Agp. 

q>-l 

For the first term we easily get that there exists a function ip £ such that 

III = j;ai{[^,(D),5,_i(£(a;e/r))]A,p}, 

with Tpq = 2^'?'(/;(2'?-). By using the Bernstein inequahty, this yields 



di{[tP,{B),Sg^i{C{LOe/r))]A,p} 
Thanks to Lemma 12.71 and (|24p , we find 



L2 



< C2'i 



[V',(D),5,_i(£(a;e/r))]A,p 



L2 



^l{[^l>giD),S,^liCiu;e/r))]A,p} < C2''||xVg||Li ||V£(a;e/r)|U^ ||A,p||i2 



< C\\xiI)\\li \\uJg/r\\L3,l \\Aqp\\L2. 



It follows that 



IIIill 1 < ^^^2"^ 9i{[V^,(D),Vi(^(W0)]A,p} ^<C\\ue/r\\L.M\ 1 

^2.1 ggf^ L B^ , 



and hence by using the embedding i?! 1 ^ -^^'^ (see Proposition 12.51 - (4)), we obtain 
(33) llnilUa.i < C\\u;e/r\\L3,i\\p\\ 



1 . 



To estimate the term II2 we do not need to detect cancellation in the structure of the commutator, 
we just write 



II2 = Y.^^'^i^[^<i(^('^s/r))Sg-ipj -Y,^l{A,{C{ue/r))n^JS,^lp}. 

q>0 q>0 

A useful remark is that thanks to the Bernstein inequalities and (j25p . we have 
(34) \\A,Cf\\L. < 2-2'?||v2£A,/|Up < 2-2'?||/||^p, > 0, p g]1,+oo[. 

This yields by using the Holder inequality and Proposition l2.5l -(3) that 

IIII2II 1 < J]2l'^||A,(£(^,/r))ViP|L2 + E2'^ll^.(^(^^^A))^«iViPllL2 



B^ 



q>0 



q>0 



< Y.2'^l\\A,C{ioe/r)\\L3{\\S,.ip\\^, + \\n,jSg.,p\\^,) 



q>0 



< ||^e/r||L3^2-«i||5g_ip| 



L6 



q>0 



q>0 k<q-2 



< 



l^eAllL^llpll 1 • 
^21 
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Hence we get from Proposition l2.5l -(4) that 

(35) I|ll2||i3,i < C||a;0/r||i3,i||/j|| i • 

^2,1 

For the term II3 we write 

II3 = 9l^[7^iJ,A,(£(a;e/r))]A,p + 5l [7^,,•, A,(£(6^e/r))]A,p 

i2>l -l<i2<0 

:= II31 + II32. 

To estimate the first term we first use the Bernstein inequahty to get 

||AfcIl3i||i2 < 2^= ^ ||[7^,„A,(£(c^e/r))]A,p||^,. 

q>k-4: 

Next, to estimate the terms inside the sum we do not need to use the structure of the commutator. 
By using again the Holder inequahty, (|34p and the Bernstein inequahty, we obtain 

\\[nij,A,{C{uje/r))]A,p\\^, < || A,(£(aje/r))||i3 ||A,p||^6 + ||A,(£(^eA))||L3 A^pH^e 

< 2-''\\ug/r\\L3\\Agp\\L2. 

It fohows by using again Proposition l2.5l -(4) that 

||II3i||l3,i < Plgill 1 <\\i^e/r\\LB E 2t(^-«)2''l||A,p||i3 <||^e/r||L3||H| 1. 

For the estimate of the low frequencies term H32 we need to use more deeply the structure of the 
commutator. We first write 

1132= Yl [dinij,Ag{C{uje/r))]\p- ^ 9l£A,(a;e/r)7^,J•A,p. 

-l<g<0 -i<g<o 

The last term of the above identity is estimated as follows by using again Proposition 12.51 -(1) and 
(3) and ([MD 

II Y diCA.iioeMTZijAgpW 1 < Yl \\diCAg{uje/r)nijAgp\\L, 

-l<q<0 "^2,1 -l<q<0 



< \\diC{ujg/r)\\L^\\p\ 
^ ll^eAllLS.illpll 1 - 



L2 



-°2,i 

To estimate the first term of H32 we write for every — 1 < g < thanks to Lemma 12.71 that 
[din.,j,Ag{C{ujg/r))]Aqp 5 < \\xh\\ w\\VC{uJe/r)\\Loo\\Agp\\L2. 

L2 1^ a 

where h{^,) = ii^^x{i) ^^id x ^ ^(K^)- Using Mikhlin-Hormander Theorem we have 

\h{x)\ < C7(l + |x|)-^ VxGMI 
This gives in particular xh G . Therefore we get by using again (f24l) that 

II Y [diTlij,Ag{C{u:e/r))\Aqp\\ 1 <\\^C{u:e/r)\\L^\\p\\L'^<\\i^elr\\L3A\p\\L^- 

-1<Q<0 

1. 

By using the embedding ^ L'^'^ which comes from Proposition 12.51 (4). we find that 

pl32||L3,i < ||a;e/r||i3,i||/)||L2. 
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We have thus obtained that the term II3 enjoys the estimate 

Consequently, by gathering this last estimate and the estimates ([33]) . (j35]) . we finally get that 



(36) 



|II||i3,i < C||a;0/r||i3,i||p|| 1 . 

^2,1 



Estimate o/III. We also decompose the term III by using Bony's formula as follows: 

III = nil + III2 + III3, 



with 



nil = di^[nij,Sg.iid3A-\uje/r))]Ag{xip) 



q>0 



III2 = di^[nij,AgidsA-\u;g/r))]Sg-i{xip) 



q>0 

III3 = diY.[ni„Ag{dsA-\ue/r))]Ag{x,p). 

q>-l 

As we have done to handle the term IIi, we can use that there exists a function £ S{W^) such 
that 

nil = ^5i([V,(D), Vi(^3A-Ho^e/r))]A,(xip)), 

q>0 

with ipq = 2^''i/'(2'-). For every p g]1,cxd[, we first write thanks to the Bernstein inequality that 



^,{[^PgiD),S,-l{^^A-\oJe/r))]A,ix^p)} 



LP 



< 02" 



[V,(D),Vi(53A-n^e/r))]A,(xip) 



LP 



Then by using successively Lemma 12.71 and the continuity of the Riesz transform (i.e. Proposition 
ES]-(3)), we get 



[V^,(D), Vi(53A-i(L^e/r))]A,(xip) 



< 



LP 



\x%\\Li\\Sq-i{Vd-iA \oje/r))\\Lv\\Aq{x 



iPn\L° 



< 2-^x^\\Li\\Vd3A-'{ug/r))\\Lp\\A,{xip)\\L^ 

< 2''^\\u)g/r\\Lp\\Ag{xip)\\Lo^ . 

It follows that 

PllilUf ^'^\\^o/r\\Lp\\Ag{xip)\\L^ < \\ujg/r\\Lp\\xip\\^<}^^ 

q>0 

This proves that the linear operator T 

f ^Y.di{[M^), 3,-1(83 A-'f)]Agixip)] 

q>0 

is continuous from L'f into itself for every p g]1,oo[ and that 

ll^ll/:(LP) < CpllxipW^o^^^. 

Consequently, by using the interpolation result of Theorem [231 we get that T is continuous on L^'^ 
for every 1 < p < 00 and q G [l,oo]. In particular, this yields 



(37) 



nil 



< 



\ug/r\\L-i,i\\xip\\QO 

00,1 
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For the term III3, we use split it into 

III3 = ai J^[7^i,•,A,(a3A-l(c^e/r■))]Aq(xlp) + al [IZij, ^^{d^^-\uoe/r))]K^{xip) 

q>l -i<g<o 
:= III31+III32. 

Let p g]1,oo[ from the Bernstein inequahty, we have that 

IIAfclllsiliLp < \\[nij,\{dz^-\^e/r))]Kq{xip)\\^^ 

q>k-A 

and the terms inside the sum can be controlled without using the structure of the commutator. We 
just write 

\\[n,j,^,{d:i^'\ue/r))]Kq{xip)\\^^ < ||A,(93A-i(u;e/r))||Lp||A,(xip)||^^ 

+ \\Ag{d3A-\u}e/r))\\Lp\\nijAg{xip)\\L^ 

< 2"'^||a;e/r||Lp||Ag(xip)||Loo. 

Note that we have used the Bernstein inequality, the continuity of the Riesz transform on and 
the fact that the support of the Fourier transform of Ag does not contains zero which gives that 
the operator TZijAg also acts continuouly on (since it can be written as the convolution with 
an fonction). It follows that for every p +00 [, we have 

IIIII31IILP < \\uj0/r\\LP ^ 2'''''\\Ag{xip)\\L^ <\\uje/r\\Lp\\xip\\^i>^^. 

k>-lq>k-4: 

By using again the interpolation result of Theorem 12.41 this yields 

l|ni3i||L3,i < \\uje/r\\j:^3,i\\xip\\BQ^^^. 

We can also estimate the term III32 without using the structure of the commutator. By using the 
continuity of the Riesz transform on and (llSp , we obtain 

\\[nij,Ag{d3A~\uje/r))]Ag{xip)\\^, < \\Ag{d3A-\uje/r))\\L^\\Ag{xip)\\^, 

+ \\Ag{d3A-\u;e/r))\\L^\\'R.ijAg{xip)\\L2 

< \\d3A-\u;e/r)\\Lo^\\Ag{xip)\\L2 
^ \\^e/r\\LB.A\xip\\L2. 

Therefore we get 

PII32II 1 < \\u}0/r\\LS,i\\xip\\L2. 

Consequently we obtain 

(38) ll™3llsii ~ ll^eA||L3,i||xip||50^^^ni2. 

Let us now turn to the estimate of the term III2. We write 

III2 = J^[7^,,•,A,(al3A-l(L^e/r•)]5,_l(xlp)+ [7e,,-,A,(93A-nt^eA)]9iVi(^i/') 

= III21+III22. 
We have by definition of the paraproducts that 

III21 = nijiT^^pUviiu^e/r)) - ^,^^^.(^^p)7^l3(a;0/r). 
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Thanks to Proposition 12.51 we get that 

I|III2i||l3,i < \\TZiz{^^e/r)\\Li,i{\\xip\\L°o + \\TZu{xip)\\l^) 

^ W^e/rWi-i^ {\\xip\\lo^ + \\Rn{xip)\\L^) 

^ ll^eAllL3.i|ki/5|lsj^ jnL2- 

Note that the norm in the right hand-side comes from the low frequency term in the Littlewood- 
Paley decomposition: we have 

(39) ||7^l3A_l(xl/))||Loo < ||7^l3A_l(xlp)||i2 < ||xip||i2 

thanks to the Berntein inequahty and the continuity of the Riesz transform. 

For the estimate of III22, we shall use that thanks to the Bernstein inequality, we have for every / 

that, 

IIA.aaA-VlUp < 2-'?||/|Up, > 0, p g]1,+oo[. 

This yields 

||ni22||LP < \\^e/r\\Lv^'i~'^{\\dlSq^i{xip)\\L'^ + \\diSq^ilZij{xip)\\L^) 

< \\uje/r\\LpYl Yl 2P-'?(||Ap(xi/9)||L- + ||Ap7^iJ(xlp)||L^) 

q>0 cj-2>p>-l 
\\^9/r\\Lp{\\xip\\B0 , +\\T^ij{xip)\\BO ) 

^ 00,1 CXD,1^ 

||^0/?-||Lp|klPllB^;^_,nL2 

by using again p9]) . Consequently, by interpolation, we also find 

I|III22||l3,i < ||a;e/r||i3,i||xip||B0 nL2. 

We have thus shown that 

(40) l|ni2||L3.i < ll^eAllL3,i|kl/3|li30^^^nL2- 

Gathering ([37]), ([Ml) and (00]) , we obtain 

(41) l|ni||L3,i < ||t^eAllL3,i||3;i/)||50^^^nL2- 
Finally we obtain 

||5i{['^ij>^^^]p}|L3,i < ll^e/?'||L3,i(||xi/9||B0^^nL2 + IIpII 1 ) 

°°' ^2,1 

thanks to ()4ip . (j36p . (j32p and (|30p . In the same way, we also obtain the estimate 
||92{[7^^J•,^;^]p}||^3,l < ||^^0/r||L3,i (HxspUso^ ^ni^ + IIpII 1 )• 

In view of (j28p . it remains to estimate the term c?3([7^jj, which has a different structure. • 
Estimate 0/ ^3 ([T^jj, w^]/?) . Since we can write that 

Av' = -(curl a;)3 = -{drUJe + u^e.) = -{rdr{^) + 2^^) = -x^ ■ Vh(.-) - 2-, 



< 

< 



— ) + 2 — ) = -Xh ■ \/h{ — ) - 2— 

we obtain that 

^^))_2A-i(^^ 



^;3(x) = -A-i(x,-V,(^))-2A- 
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and hence by using Lemma 12.101 that 

2 

-v^{x) = Xh-VhA-\ue/r)-2Y^A-^TZii{u}e/r) + 2A-\uje/r) 

i=l 

(42) = Xh-VhA-\coe/r) + 2A-^n3s{uje/r). 

Thus, we have a decomposition of the commutator under the form 

2 

-d3([nij,v'^]pj = Y,^3(dkA'\uje/r)[n,j,Xk]p^ +2d3([n,j,A-^ns3{uJe/r)]i 



k=l 
2 

+ ^d3(^[n^j,dkA-\oje/r)] (xkp) 

k=l 

T + IT + III. 



To estimate the first term I, we use Lemma l2.10l -(2) to obtain that 

ds[dkA-Hy)[n,j,Xk]p) = d:,(duA-\^)c%p) =n^kC^) cy + dkA-\^)d,C%p. 



r 

It follows that 



|I||l3,i < (ll4p||L-||^3feK/r)||L3.i + \\dkA-\,oe/r)\\L^\\d^C%p\\LZ,. 

k=l 

^ IIpIIl3.i||w0AIIl3.i 



thanks to (|26p . (j27p . The estimates of the terms II and III are similar to the ones of II and III 
in (pO|) (indeed, the operator A~^7^33 = d^^A^'^ has the same properties as £ = — 25i3A^^ which 
arises in ([30]) ) consequently, we also get as in ([36|) and ([^T|) that 



|I||l3.i < ll^eAllLsllpll 1 , ||n||2,3,i < \\u}e/r\\L3,i\\pxh\\BZ 

Br, , 



Consequently, we also find that 

||53([7^^J•,^;3]p) 11^3,1 < ||w0/r||i3,i (llpx/ill^o^ ^nL^ + IIpII i )• 

This ends the proof of Theorem 13.11 □ 

3.2. Commutation between the advection operator and Aq. The last commutator estimate 
which is needed in the proof of our main result is the following. 

Proposition 3.2. Let v be an axisymmetric divergence free vector field without swirl and p a 
smooth scalar function. Then there exists C > such that for every g G N U { — 1} we have 

II [Aq,v- V]p\\^2 < C\\uje/r\\Ls.i{\\pxh\\L6 + \\p\\l^)- 

Proof. From the incompressibility of the velocity we have 

3 

(43) [Aq,vV]p = Y.di(^[Aq,v']p) =1 + 11 + 111. 

i=l 

The first and the second terms can be handled in the same way, so we shall only detail the proof 
of the estimate of the first one. Thanks to (1291). we have that 



v^{x) = xiA-^dsiug/r) + C{oje/r), with C = -27^13A"^ 
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and hence we get 



di 



Ag,xiA-^d3iu;g/r) p) + d_ 



Ag, d^^e/r) 



h+h 



The estimate of the second term in the right-hand side is again a direct consequence of Lemma 
and ()24p . Indeed, we write 

P2IIL2 ^ ||V£(u;eA)||^ooll/5||L2 ^ ||^eA||i3,ill/'llL2. 
The first term Ii in the right-hand side can be expanded under the form 

Ii = di 



Iii+I 



12- 



A„ A-^d3{ue/r)\ (xip)) + l^A-^dsicog/r) [A„ xi] p 
We start with the estimate of Ii2- By definition of Ag, we have 

xiAgp = xi239 /" ^{2'i{x-y))piy)dy 

= 23'?/ ^{2'i{x-y))yip{y)dy + 2^'i I ^{2'i{x - y)){xi - y^)p{y)dy 

= A,(xip) + 2-^23^1(2"-) 
where ^i{x) = xiip{x) G Consequently we get the expression of the commutator: 

(44) [A„xi]p = -2-^23^1 (2«-)*P- 

This yields 

I12 = -{nr3{u;e/r))2^''ipi{2'i-)*p-{A-'d3{u:e/r)}2^Hdiipi){2''-)^p. 

Therefore we get by using again the Holder inequality, the continuity of the Riesz transform, ([18 
and the Young inequality for convolutions that: 

||Ii2||l2 < \\nn{oJe/r)\\^,2^''\\M'^''-) * p\\^, + \\A-'d5{iJe/r)\\^^2^'^ 

^ l|we/r||L3||99l||^3 ||/>||i2 + ll^^e /r 11^3.1 ||(9l 11^1 ||/)||i2 



< 



ll'^e/?'||L3,i 11^11^2. 



Note that we have also used the embedding ([15 
To estimate In we use again Lemma l2.8t 



iii||L2 l|VA ^93(we/r)||i3||a;ip||i6 < ||a;0/r||i3||xip||i6 < ||tj0/r||i3,i||xip||i6. 
We have thus shown that 

(45) ||I||l2 < ||a;e/r||i3.i {\\p\\lp + Wxiph^^)- 
In the same way, we obtain that 

(46) ||II||l2 < \\uje/r\\L3,i{\\p\\L2 + \\x2p\\l6)- 
It remains to estimate the last term III. By using (j42p . we get 

-III = d3{[AgyhA-Hu;e/r)]{xhp)}+d3{VhA-\ioe/r)[Ag,Xh\p} 

+ 293{[Ag,A"i7l33(we/r)]p} 
= nil + III2 + III3. 

The estimates of the first and last terms follow again from Lemma 12.81 we write that 

||nii||i2 < C||V^A"^(we/r)||i3||x/,p||i6 < \\uj0/r\\L3\\xhp\\LS < \\u;e/r\\L3,i\\xhp\\Le 
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and that 

||ni3||2,2 < C||VA"^7^33(a;0/r)||Loo||p||^2 < \\^Zz■i{uJe/r)\\L3^\\p\\L2 < (:7||a;0/r||i3,i ||p||l2. 
Note that we have used again the estimate (jlSp . 

Finally, to estimate the second term III2 we can use the expression of the commutator [Aq,x/i] 
given by 



III2 = 2-^d3(^{VhA-Hoje/r))2^\h{2'^-)*p'^ 

= 2-'^{d3VhA-\u;e/r)) {2^'^ipf,{2i-) * p) +VhA-\ue/r) {2^'i{d3^h){2''-) ^ p) , 
with (ph{x) = —Xhfix). It follows as before that 

IIIII2IIL2 < 2-^||a;,/r||^323^||(^,,(2''-)*p||i6 + ||V;,A-i(^e/r)||Loo23'?||(a3(^0(2''-)*/'llL2 



< 



\\'^e/r\\L3A\\p\\L^- 
Gathering these estimates we also find that 

PII||2,3 < \\uJg/r\\L3^{\\p\\LP + WxhpWL^)- 

In view of (P3l) . (Il5]l . (fl6l) and the last estimate, this ends the proof of Proposition 13.21 □ 

4. A PRIORI ESTIMATES 

In this section we intend to establish the global a priori estimates needed for the proof of 
Theorem ll.il We shall first prove some basic weak estimates that can be obtained easily through 
energy type estimates. In a second step, we shall prove the control of some stronger norms such as 
||w(t)||L°° and ||Vu(t)||Lcx) . This part requires more refined analysis: we use the special structure of 
the Boussinesq model combined with the previous commutator estimates. 

4.1. Energy estimates. We start with some elementary energy estimates. 

Proposition 4.1. Let {v,p) be a smooth solution of ([1]) then 

(1) For p £]l,oo[,q £ [1, 00] and t G M+, we have 

+ 2||Vp|||2^2 < ||po|li2, ||p||l^lp.<? < C'||/?o||lp>'?- 

(2) For vq £ L'^,po £ L?' and t £ M+ we have 

\\v{t)\\L2 < \\vo\\l2 +t\\po\\L2- 

(3) For pq G we have the dispersive estimate 

\\pmL^<cM^- 
ti 

The constant C is absolute. 
Note that the axisymmetric assumption is not needed in this proposition 

Proof. (1) By taking the L^-scalar product of the second equation of ([1]) with p and integrating by 
parts, we get since v is divergence free that 

1 d 



2dt 



I2+ [ \Vp{t,x)\'^dx = 0. 
Jr3 



Integrating in time this differential inequality gives the desired result. 
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Let us now move to the estimate of the density in Lorentz spaces. First, the same argument yields 
that for every p € [1, oo], we have 

\\p{t)\\LP < WPoWlp- 

It suffices now to use the interpolation result of Theorem 12.41 

(2) We take the L-^-scalar product of the velocity equation with v and we integrate by parts 



and this implies that 

Thus, integrating in time gives 



L2. 



|l2 < \\vo\\l2 + / \\p{T)\\L2dT. 

Jo 

Since ||p(t)||j;^2 < ||po||l2, we infer 

||t'(i)llL2 < \\vo\\l^ +*II/5o||l2. 

(3) The estimate is a direct consequence of Lemma fA.ll The proof of the proposition is now achieved. 

□ 

4.2. Estimates of the moments of p.. We have seen in subsection 13.11 and subsection 13.21 that 
the estimates of the commutators involve some moments of the density. Thus we aim in this 
paragraph at giving suitable estimates for the moments that will be needed later when we shall 
perform our diagonalization of the Boussinesq system. Two types of estimates are discussed: the 
energy estimates of the horizontal moments \xh\^p, with k = 1,2 and some dispersive estimates. 
More precisely we prove the following. 

Proposition 4.2. Let v be a vector field with zero divergence and satisfying the energy estimate of 
Proposition\4. 1\ Let p be a solution of the transport- diffusion equation 



dtp + v -Vp- Ap = 0, p{0,x) = po. 
Then we have the following estimates. 

(1) For po £ o,f^d x^Po G L^ , there exists Cq > such that for every t > 

(2) For pq £ L'^ n L™", m > 6 and x^po G L'^ , there exists Cq > such that for every t > 

\\xhpit)\\L^ < Coit^+t'i). 

(3) For pq G L^ and \xh\'^pQ G L^, there exists Cq > such that for every t > 

llk/iP/0||Lj°°L2 + \\\xhfp\\L2H^ ^ '^0(1 + t^)- 

(4) For po G -L^ n L^ and \xh\'^Po G L^, there exists Cq > such that for every t > 

\\\Xh\^p{t)\\L<i < Co{tf +t-'2). 

Remark 4.3. Note that when pQ G L^ and |x/ip/9o G L^ then automatically the moment of order 
one belongs to L^, that is XhPo G L^. This is an easy consequence of the Holder inequality 

1 i_ 

Wxhph^ < \\p\\l2\\\Xh\'^P\\l2- 
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Proof. (1) Setting / = XhP, we can easily check that / solves the equation 

(47) dtf + v-Vf-Af = v'^p - 2VhP 

with the notations = {v^ ,v^) and V/^ = (5i,52). Now, taking the L^-scalar product 
with /, integrating by parts and using the Holder inequality 

^:|ll/wili^ + liv/(t)iii. = / v^pfdx-2j v.pfdx 

< \ML4p\\L4f\\L^ + np\\L4^f\\L^- 

By using the Sobolev embedding ^ combined with the Young inequality, we obtain 

|ll/(i)lli2 + iiv/(t)iii.<ibiii.iHii3 + iiHli.. 

Since the Gagliardo-Nirenberg inequality gives that 

WpWls < II/'IIl2||Vp||j,2, 

we infer 

(48) + liwwili^ ^ Mh\\p\\L4^p\\L- + \\p\\h- 

Integrating in time and using the energy estimate of Proposition I4.H -(1). we thus obtain 

Wfmh + l|V/(t)|li2^2 < \\M\I2 + \\v\\Io.l4po\\l4^p\\lIl^ + WPoWl^t 

< WfoWh+Coil + tY^ + \\po\\l2t 

< Co(l + ti). 

(2) We shall apply Lemma lA.ll to (I47p with F = pti and G = Vip. First, we observe that we have 
obviously from the Holder inequality combined with Proposition 14 . 11 - ( 1 ) that for m>2 

W^Wl^l^ ^ II^^IIl-l^IIpIIl-l- 

< Co(l + t) 

and 

ll-^llLfL6 < IIpoIIlS- 

Consequently, we get from Lemma lA.ll and Proposition 14.11 that for m > 6 and for t > 0, 

||/(t)||L- < C(l+t-t)||/o||i2+Co(l + t3-2l7) + (l+t3)||po||^6 

< Co(i-3+t3). 

(3) The second moment g = \xh\^p solves the following equation 

dtg + v-Vg- Ag = 2v^{xhp) - 2VhP - MWh{xhp) 
= Vhf -2Vhp- ^divhf. 
By using again an energy estimate, we find that 
d_ 
It' 

Thus, by integrating in time and using the energy estimates for p, v and / we get 

II^WIIl^ + \\Vg\\L2L2 < \\go\\L2 + \\v\\LrL4f\\l^L4'^f\\l2L2t-^ + llpolli^i^ + II/IIl-l^*^ 

< Co{l + tl), 

where Co is a constant depending on the quantities |||x/j|'^'po||l2 for A; = 0, 1, 2. 
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znumh + wvgmh < \\vmi2\\f{t)h2\\vf{t)h2 + mmh + wmwh. 



(4) By setting gi{t,x) = tg{t,x), we have that 

dtgi+v -Vgi- Agi = g + 2v^{txhp) - 2tV hp - ^^^^h{txhp). 

Multiplying this equation by l^il^^i, integrating by parts and the obvious inequality \xhp\ < 
|/)| 2 l^f]^! 2 J we thus get 

^^ll5i|li6+5 / iV^il^l^il'd^ < NlL6|bi|li6+t^ / \v\\p\^g,\'^ dx 

+ t! \p\\VgMUx + t"2 [ \p\"2\g,\l\Vgi\dx. 

It follows from Holder inequality that 

[ \p\h\g,\^Vgl\dx<t'2\\v\\L2\\g,\\fM\Ks 

Consequently 

ItMWI^^ + ^I \^9i?\9i\'dx < NlL«lbi|li6+t^|HU2||ffi||^8||HlJ^ 
D at 

+ m\L49i\\h+t'^\\p\\U9i\\h)( [ \Vg,\'\g,\'dx)\ 
Now we can use the Young inequality combined with the following Sobolev inequality 

llffi|li'i8<l|V(<7?)||i2=9 / \Vg,\'\g,\Ux 

to obtain that 

yM\\% + 49i\\U + c [ \Vg,\'\g,\Ux < \\9\\L49i\\h + t%v\\fM' ^ 
dt J ^3 

+ t4p\\l49i\\\^+t\\p\\L49i\\h- 

By using Proposition 14.11 we deduce that 

;^ll5i|li6 + c|bi||ii8+c/ \Vg,\'\g,\'dx < \\g\\^4g,\\l, + Cot%l + t'') 
dt Jr3 

+ t4po\\U9l\\h+t\\po\\L49l\\le. 

Next, by using again the Young inequality, we infer 

|||5l(t)|lia < Co(i6+P)+Co(t+||9(t)||LOlblWlli6. 

By integrating in time this differential inequality, we obtain that 

\\9i{t)\\% < Co{t' + t'') + Co [\r + \\giT)Ue)\\gi{T)\\%d^ 

Jo 

Therefore we get from Proposition l4.2l -(3) combined with the Sobolev embedding C that 

7 19 

||9i(t)llL6 < Co{tf> +tir) + Co\\g\\LiL6 

7 19 1 

< Co{t-^ +t-) + t2\\Vg\\L2L^ 

< Co{tl +tf) + Co{l + t^)t^. 

Therefore, we obtain that 

\\\xH\'pmL<'<Co{t'i 
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This ends the proof of Proposition 14.2 



□ 



4.3. Strong estimates. As in the study of the axisymmetric Euler equation, the main impor- 
tant quantity that one should estimate in order to get the global existence of smooth solutions is 
||^(t)||i3,i. Indeed, this will enable us to bound stronger norms such as ||a;(t)||j;,oo and ||Vt'(t)||Loo 
which are the significant quantities to propagate higher regularities. 

4.3.1. Estimate of \\^{t)\\Loo. First, we will introduce the following notation: we denote by any 
function of the form 
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^kit) = Co exp(.^..exp (Cot^)...), 

k times 

where Cq depends on the involved norms of the initial data and its value may vary from line to 
line up to some absolute constants. We will make an intensive use (without mentionning it) of the 
following trivial facts 

[ <^k{T)dT < <^k{t) and exp( [ ^k{r)dT) < <^k+i{t). 
Jo Jo 

We first establish the following result. 

Proposition 4.4. Let vq G be an axisymmetric vector field such that ^ G L^'^ and po G L^n-L™, 
for m > 6, axisymmetric and such that \xh\'^Po S L^- Then, we have for every t G M_|_ 

ll7WlLM + llyWlLoc<^2(i), 

where Cq is a constant depending on the norms of the initial data. 

Proof. Recall that the equation of the scalar component of the vorticity uj = ujocq is given by 



(49) 



dtUjQ + V ■ VcOg = — UJQ — drp. 

r 



It follows that the evolution of the quantity — is governed by the equation 



(50) 



{dt + v V] 



drP 



By applying the operator ^ to the equation of the density in ([TJ, we obtain that 

(dt + v V) i-drA-'p) - ^ = - \-drA-\v ■ Vl p. 

By setting T := + p, vfe infer 



{dt + v-V)T 



Observe that the incompressibility of the velocity field allows us to get that for every p G [1, cxd] 



|r(t)||LP < ||ro||LP + 







p 



dr. 



LP 



Therefore we get by the interpolation result of Theorem 12.41 that for 1 < p < oo and q & [l,oo] 



LP-1 



^ l|rol|LP'9 + 



Lr 



P{r) 



LP'i 



dr. 



In particular, we have 



L3.1 < l|ro||L3,i + 



/' 

Jo 



1 



drA-\vV p{t) 



L3,l 



dr. 
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Applying Theorem 13.11 we find 



|r(t)||L3.i < ||ro||i3,i+ / ||(weA)(r)||i3,i(||x,,p(r)||B0^^nL2 + 11/^(^)11 1 )dr. 



2,1 



Moreover, thanks to Proposition 12.91 and Proposition 14.11 we have 

||(^,/r)(t)||i3a < ||r(t)||i3a + ||^5,A-V(i)L3.i <l|r(t)||L3,i+C||/5o||i3,i. 
The combination of these last estimates yield 

ll(^eA)(t)||L3,i <C{\\ujo/r\\L3,i + \\po\yA)+ [ \\{^e/r){T)\\L3,i{\\xhp{T)\\j,o 0^2 + ||p(r)|| i )dr. 

Jo °°- ^2,1 

Thus we get by the Gronwall inequality that 

(51) \\{u}e/r){t)\\L3,i < C(||woAllL3a + ||/9o||L3,i)exp(C||x/,p||j^i(Bo^^nL2) + ^11^11 1 ). 

The term ||p|| 1 will be controlled only by energy estimates. Indeed, the interpolation estimate 

■^(-82,1 



pII 1 < IIPIII2IIVHI 



L2 



combined with Proposition 14.11 and the Holder inequality give 



1 , 1 



IIpIL^I < \\p\\l^L^t'^\\'^P\\l2L2<t~^\\po\\L^- 

To control the term ||2;/ip||iii2 in the right hand side of (|5ip . we can use Proposition 14.2 

9 

\\XhP\\L}L'^ < t\\Xhp\\Lr'L2 <Co{l + ti). 



Consequently we obtain in view of (151)1 

(52) ||^(t)|L3.<Coe^°*^e^"^^^"-*^-^.. 

Now it remains to estimate the right term of (I52p inside the exponential. Let us first sketch the 
strategy of our approach. We will introduce an integer N{t) € N that will be chosen in an optimal 
way in the end and we will split in frequency the involved quantity into two parts: low frequencies 
corresponding to g < N{t) and high frequencies associated to g > N{t). To estimate the low 
fequencies we use the dispersive result of Proposition l4.2l -(2). The estimate of high frequencies is 
based on a smoothing effect. 

By using Proposition l4.2l -(2) and the Bernstein inequality, we find that 

W^hphlBO^^ = [ J2 \\^q{Xhp){r)\\L°-dT + j ^ || Ag (xfep) (t) ||loo 
°°' q<N{T) q>N{T) 

(53) < Co [\t-^ +T-i)N{T)dT + C f V 2«i||A,(x/,p)(T)|U2dr. 

q>N(r) 

Now we intend to estimate the last sum in the above inequality. For this purpose we localize in 
frequency the equation for f = XhP which is 



dtf + v-Vf-Af = v^p - 2VhP ■■= F 



By Setting /„ := A„/, we infer 



dtfq +V-Vfg- Afg = -[Ag, V ' V] f + Fg 
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From an L? energy estimate, we obtain that 

\jt\W)\\l- - Ij^f'^'^f'^^'' ^ ||/,||l.(||[A„^ • V]fh2 + \\F,U2] 
Since the Bessel identity yields 

c2'"^\\U\l2 < - [ {M<i)hdx, 



it follows that 



dt 

Therefore we obtain by integration in time that 

ct 



IIAWIIl^ <e-^*2'l/,(0)||^2+ / e-^(*-^)2'^(||[A„T;.V]/||^. + ||F,||^.)dr. 

JO 

To estimate the commutator in the right hand side, we can use Proposition 13.21 and Proposition 14.2 
ll[A„^;-V]/(r)||i. < C||(a;,/r)(r)||^3,i(|||x,|V(r)||i6 + ||x,/,(r)||^.) 



< Co||(a;e/r)(r)||i3,i(Tf +T-^). 
Hence we get 



Jo 

(54) + Co [\-<'-^^^"''\\{ue/r){T)hs^{Tf +T--2)dT. 



13 1 

Let us set /C(r) = t~ + , then ([5^) and convolution inequalities yield 

/* 2''i||A,(x,p)(T)b.dr< J] 2-5'^(||/,(0)|U. + ||F,||^i^.) 

q>N(T) q>-l 

+Cof Yl re-(---')2^'/c(r')||(u;e/r)(r')||i3,idr' 

q>N{r) 

q>N(T) ■'^ 

Moreover, from Proposition 14.11 and Lemma lA. 11 we also have 

\\F\\lIl^ < t'2\\Vp\\L2L^ + \\v\\L^L4p\\LjL^ 

< C\\poh2th +Co{l + t)\\po\\L2(^j^ r-idr + i) 
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Inserting these estimates into (|53|) yields 

\\xhp\\L]B^ , < Co{l + t^) + Co ['{t-* +T-'i)N{T)dT 

oo,i 

/■* 1 -i 

< Co{l + t^) + Co (t4 +T-4)iV(r)(iT 
Jo 

(55) + Co^Ve/HlL-L3.i(T'#2-5^M+ J] 2"! J\-<--'^^"'{T'}-"^dT')dr. 



g>N{T) 

By a change of variables we get 



q>N{T) ■'^ 



g>Af(T) 



I(r)+II(r). 



with 



Ai(r) = |g > iV(r) and r22'' > l| and A2(r) = |g > iV(r) and r22'? < l| . 



To estimate the first term we use the following inequality which can, be proven by integration by 
parts: there exists C > such that for every x >\ 

r -1 _i 
/ y 2e'^ydy < Cx 2e'^^. 

Jo 

It follows that 

q>N{T) 

To estimate the second term, we observe that the integral is bounded by a fixed number and hence, 
we find that 

II(r) < 2^'^<2-^^M Y 2''<2-5^M(l + r-^). 

i}gA2(r) 22'J<r-l 

Gathering these estimates, we obtain 

Y 2^^ re-^(^-^')2''{r'}-^(ir'<2-5^M(l + r-5). 

q>N(T) ''^ 

By plugging this estimate into ([55]) , we get 



\\xhp\\ 



,1^1^^ < Co{l + t^) + Co (Tf+T-t)(iV(r) + ||a;e/r 11^00^3,1 2-i^M) dr. 
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We choose such that 

N{t) = 2[log2 (2 + \\ue/r\\L^L3.i) 
and then we find 

(56) \\xhP\\L}BO , < Co{l + tf) + Co (rf + t'^) log (2 + \\ujg/r\\L^L3.i)dT. 

Putting together ([52]) . ([56]) and Proposition l4.2l -(2). we find that 
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log(2 + ||we/r||rooi3,i) < Co(l + t~) + Co / {t- + T~^)log{2 + \\uje/r\\LooL3,i)dT.. 

Jo 

From the Gronwall inequahty, we infer 

log (2 + ||a;e/r||ioo^3,i) < Co{l + tir)e^o(i^+t^) < .^^^(i). 



Therefore we get by using again (j56j) that 

II 



r 



Since ^ = ^ee, impHes that 



(57) ||^(t)||^3,i<^2(t). 



Finally, thanks to p9|) . we obtain 

llyWIlLoo <C||^(t)||z.3,i <f2(0- 

This ends the proof of Proposition I4.4[ □ 
4.3.2. Estimate of \\uj{t)\\L°° ■ Our purpose now is to bound the vorticity. 

Proposition 4.5. Let vq G be an axisymmetric divergence free vector field without swirl such 
that ujQ G L°°, ^ G L'^'^. Let po axisymmetric scalar function, belonging to L^ riL^,m > 6 and 
such that \xh\'^po G L^. Then we have for every t G 

Mt)\\L^ + \\yp\\LlL^<Mt)- 
Proof. From the maximum principle for the equation we obtain that 



||w(t)||Loo < ||woI|l°° + / Wv"" /r{T)\\L^\\uj{T)\\L^dT + / \\Vp{T)\\L^dT. 

Jo Jo 
By combining Proposition 14.41 and the Gronwall inequality, this yields 

Now we claim that, 

(58) \Np\\LlL^<Co(l + t^ + \\u{T)\\LoodTy 

Let us first finish the proof by using this estimate. We deduce that 

II- < $3(t)(l + ^ \\uj{T)\\L^dTy 



and thanks to the Gronwall inequality that 



Mt)\\L^ < Mt)- 
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ndX. 



LP + 



This gives in turn 

||V/5||^lioo < $4(t). 

Let us now come back to the proof of ([58]) . For g G N we set pq := A^p, then 

dtPq + V ■ Vpg - Apq = -[Ag, V ■ V]p 
Let p >2 then multiplying this equation by \Pq\^~'^Pq and using Holder inequality 

-YM^Wl.- [ {Apq)\pqr'pqdx<\\pq\\l-'\\[Aq,V-V]p\\^^. 

Now we use the generalized Bernstein inequality, see [23] . 

-2'''\\Pq\\l,<- [ (Ap,)|p,rv 

Hence we get 

d 

di' 
This gives 

(59) IIp.WIIlp <e-^-*2^lA,po||Lp+ T e''^''^''^*-^) || [A„ t; • V]/>||^,dr. 

JO 

Integrating in time implies that 

MlIlp < 2"''IIpo||lp + 2-2^11 [Ag, • V]p\\^,^^ 
According to Proposition 2.3 [17J and Proposition 14. II we have 

\\[Aq,vV]p\\^^ < C\\p\\Lp{iq + l)\\uj\\L^ + \\VA^iv\\l^) 

< C\\po\\Lp{{q + 1)||w||l°° + \\v\\l^) 

< C||po||Lp((g + l)||w||L- +Co(l + t)). 

It follows that 

WPghjLP < Co{l + t^)2-^'^ + C\\po\\Lv{q + 1)2-2511^^11^1^^. 
By using the Bernstein inequality, we find for p > 3 that 

llVpll^i^oo < Ct||/>o||L2 + Cj;2^(^+|)||p,||^i^, 

< Co(l + i^) Yl 2^^"'^'^ + C^ollo^ll^jz^o. 5] 2«(-i+|)(g + 1) 

< Co{l+t^) + Co\\io\\LiLoo. 

This ends the proof of the desired inequality. □ 

4.3.3. Lipschitz bound of the velocity. We shall now deal with the global propagation of the sub- 
critical Sobolev regularities. This is basically related to the control of the Lipschitz norm of the 
velocity. 

Proposition 4.6. Let | < s < 3 and (fo,po) G H'^xH'^^^ and {v,p) be a solution of the Boussinesq 
system ^ . Then we have for every i > 

II II I II II I II II ^ n ^^ +\ '^ll^^lliiioo 

If in addition pQ G L™ with m > 6 and \xh\'^po G L^, then we get for every t > 

''^lljooffs + 11/^112^00^:^8-2 "I" ll/Olljljl^s ^ *1'6(^)- 
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Remark 4.7. We point out that we can extend the results of Proposition 14.61 to higher regularities 
s > 3 but for the sake of simplicity we restrict ourselves here to the case of s < 3. 

Proof. We localize in frequency the equation of the velocity. For g G N U {— 1} we set Vg := AgV 
and pg := Agp. 

dtVg + V ■ VVg + VVTg = ^gC^ " [Ag, V ' V]V . 

Thus taking the L^-scalar product with Vg and using the incompressibility of v and Vg we get 
Integrating in time we obtain 

\\Vgit)\\L2 < \\Vg{0)\\L2 + ||Pg||ili2 + \\[Ag,V V]t;||ilj^2. 

Thus we get 

t t 

We will use the commutator estimate, see for instance Lemma B.5 of [9], 







Putting together these estimates and using Gronwall inequality yield 

(60) IHl'th- ^ {\\vo\\h- + IIpIIzihO^^*^"^""'''''"- 



Using the estimate (j59|) we get for g G N 

WpghrL^ + '^^IPghjL^ < C\\pg{0)\\L2 + \\[Ag,V VjphiL^. 

Therefore we find 

\\p\\lo.Hs-2 + WpWlIh^ < \\^-IP\\lIL^ + II/>0||h-2 + ||(2''(^-2)||[^^^^ . V]p\\LlL^)g\\i2 

< Ct\\po\\L2 + \\po\\h^-2 + II (2^(^-2) II [A„ • V]/5||^ii2)g||,2. 
Since — l<s — 2<1 then we have the estimate, see [9], 

II (2^(^-2) II [A,, • V]HLilO.II^^ < C ^ W IIl^ ||p(r) ||^,.-2 dr. 

Consequently, 

||/9||;^ooj^,_2 + IIpII;^!^, <Co(l + t) + C / ||Vt;(r)||L^||p(r)||^.-2(ir. 

* * JO 

By Gronwall inequality 

(61) \\p\\z^H^~2 + IIpIIzj^. <Co{1 + t)e'""'""4^°° . 

Combining this estimate with (|60p gives the desired estimates. 

Now to get a global bound for Lipschitz norm of the velocity we use the classical logarithmic 
estimate: for s > | 

I|Vi^||loo < ||t;||L2 + ||u;||l°° log(e + ||w||//0- 
Combining this estimate with the first result of Proposition 14.61 and Proposition 14.51 

||Vv||loo < $4(t)(^l + y \\Vv{T)\\Lo^dT 

It follows from Gronwall inequality that 

||Vt;(t)||Loc < $5(t). 
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Plugging this estimate into (j60j) and (|6T]) gives 

hWi^Hs + \\p\\l^h^-2 + IIpIIli//. < ^6(t)- 

This ends the proof of the proposition. □ 

5. Proof of the main result 

The proof of the existence part of Theorem 11.11 can be done in a classical way by smoothing out 
the initial data as follows 

V0,n = SnVo = 2^^x{2'' ■) * Vo, Po,n = S^PO = 2^^x{2''-) * PO 

where Sn is the cut-off in frequency defined in the preliminaries. Since x is radial then the functions 
VQ^n and /Oo,n remain axisymmetric. Moreover this family is uniformly bounded in the space of 
initial data: this is obvious in Sobolev and Lebesgue spaces but it remains to check the uniform 
boundedness of the horizontal moment of the density. We will show that 

(62) sup|||xfep/>„,o||L2 < C(||po||l2 + |||x/jP/)o||l2). 

neN 

For this purpose we write 

\xh\'^\Pn,oi^)\ = / x(2"(x - y))p(y)(iy 

< 2/ \xh-yhf\x\{2"{x-y))\p\{y)dy + 2 I \x\{2''{x - y))\yh\^p{y)dy 

< 22-2"(23"xi(2"-) * p)ix) + 2(1x1(2"-) * i\yh\'p))ix) 
with Xiix) = \xh\'^\x{x)- From convolution laws we get 

llk/i|^l/On,o||L2 < C2~'^"'\\p\\l2 + C\\\xh\'^po\\L2. 

This achieves the proof of (|62p . Now, by using standard arguments based on the a priori estimates 
described in Proposition 14. 6t Proposition 14.51 and Proposition 14.21 we can construct a unique global 
solution {vn, Pn) in the following space 

Vnee{W+-H')r\L^{W+;W^^^) and pn e C{R+- H'-^) r\ L^{R+-W^^^). 

The control is uniform with respect to the parameter n. Therefore we can prove the strong con- 
vergence of a subsequence of {vn, Pn)n&^ to some {v,p) belonging to the same space and satisfying 
the initial value problem. It remains to prove the uniqueness problem. This gives the existence of 
a solution. 

The uniqueness will be proven in the following space 

{v,p) ex ■= (c(M+;L2)nL^(M+;VF^'°°)'^ 



Let {v^, p^) e X,l < i < 2he two solutions of the system ([T]) with the same initial data {vq, 9q) and 
denote 5v = v"^ — ,5p = p^ — p^ . Then 

{dt5v + • V5v + vn = -5v ■ Vv^ + Jpe^ 
dt5p + ^2 . V5p -A6p = -6v Vpi 
divTj* = 0. 

Taking the L^-scalar product of the first equation with 5v and integrating by parts gives 
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Consequently, 

^\\Sv{t)\\L2 < \\Vv^\\l^\\Sv\\l2 + \\Sp\\l2. 
Using the Gronwall inequality yields 

By the same computations we get 

mmL^ < \\SP0\\l^+ f\\Vp\T)\\L^\\Sv{T)U2dT. 

Jo 

It suffices now to put together these estimates and to use the Gronwall inequality. 
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Appendix A. Nash-De Giorgi estimates for convection-diffusion equations 
Lemma A.l. Consider the equation 

(64) dtf + u-Vf-Af = V-F + G, t>0, x G M^, /(O, x) = /o(x). 
Consider p, q, pi, qi,£ [l,+oo],r G [2, +00] with 

2 3 2 3 

- + -<1, — + — <2. 

p q Pi qi 

There exists C > such that for every smooth divergence free vector field u, for every F G L^L'^ 
and for every /o G L"^ , the solution of (jMj) satisfies the estimate: for every t g]0,T], 

II/WIIl- < c(i + 4)ll/ol|L^ + c(i + ^/r'"^''''^)||F||i.^, 

(65) + c(i + Vt'"^^+^^)||G||^..^,,. 

Proof. Since the equation is linear, we can study separately the three problems 



(66) 



Vf = V-F, f Vf = G, S Vf = 0, 

/(0,x) = ' 1/(0,x) = 0' 1 /(0,x) =/o(x) 



where we have set Vf = dtf + u ■ V/ — A/. 

Let us start with the first problem in (166p . We shall prove that there exists C > such that for 
every F with Hi^HiP^g < 1, we have the estimate 

(67) \\f\\LrL^<C. 

Once this estimate, is proven, the estimate involving F in (I65p will just follow by a scaling argument. 
The first step is to use the standard L'^ a priori estimate (obtained by multiplying by Vf by 
|/|''~"'^sign /). Since u is divergence free, we have that 

^J-\\fmh) + iq-i) [ NfWr'dx<{q-i) [ \F\\vf\\fr'dx. 

dt\q / Jm.2 J^s 

From the Young inequality and the Holder inequality (note that since 2/p-\-3/q < 1, we necessarily 
have that q > 3 and p > 2) this yields 



j^{hfm%)<iQ-mF\\h\\f\\i-.' 

and hence by integration in time, we obtain since q > 3 that 

(68) II/IIl^l. < \\Fit)\\hdty < < C, 

where Cq depends only on q. To improve this estimate that is to go from the above L'^ estimate to 
an L°° estimate, we shall follow the De-Giorgi, Nash iteration argument. For M > to be chosen, 
let us take a positive increasing sequence {Mk)k>o such that Mk < M and converges towards 
M. A good choice is for example 

(69) Mfc = M(l-^). 
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We shall use the standard notation x+ = max(x,0). Since u is divergence free, we obtain the level 
set energy estimate 

|(i||(/-Mfc)+(t)||i.) + ||V(/-M,)+||i. < j^^^^^\F\\Vf\dx 

< (I \F\')'\\V{f-Mk)+\\L2 

^ Jf>Mk ' 

where the last inequality comes from Cauchy-Schwarz. By using the Young inequality, we thus 
obtain 

(70) t/fc < / / \Ff dxdt 

Jo Jf>Mk 

where 

Uk = IK/ - Mk)+\\l^L'^ + ||V(/ - Mfc)+||2,^,. 

The main idea is to prove that the right-hand side of (|7U|) can estimated by a power of C/fc-i strictly 
larger than 1. By using the Holder inequality, we first get that 

(71) [ ll^(*)llL"^fc(i)'"'^*< ll^llifL.(^'"^'tW^'"'^^^^^*)'~'- 
where mk{t) = f{t) > Mk}\. To estimate mk{t), we note that if f{t,x) > then 

f{t,x)-Mk-i>Mk-Mk-i>0 



and thus we have 

(72) 1/(M)>M, < ^^^(/(t,x) - Mfc„i)+. 
This yields 

(73) mfc(t) < ^ \\{f{t) - M,^,U\\f^ 

for every m > 1. We shall choose m carefully below. By plugging this last estimate in (I7ip . we get 

(.4) < ii^iii„„((i±i)!)"'"-f (yf' um-M.^.uC-'"^T' 

Now let us notice that if a > 1 and /3 G [2,6] are such that § + ^ ^ | then we have 

\\{f{t)-Mk.i)+\\l^L,s<Uk^,. 

Indeed the control of U^-i gives a control of the L^L^ and the L\H^ norm. By Sobolev embedding 
this gives a control of the L\L^ norm and then the inequality follows by standard interpolation in 
Lebesgue spaces. Consequently, to achieve our program, we need to choose m G [2, 6] such that 

(75) ^(1_2)>2, 2 + 3 > ^ 



rail - ^ 771-2" 



The first constraint can be satisfied as soon as2/(l — 2/g)<6 which is equivalent to g > 3 while 
the second constraint can be satisfied as soon as 

2(1 - -) + 3(1 - -) > m(l - -) > - • 2 = 3 
p q q 2 

which is equivalent to 

2 3 

- + - < 1. 

p q 
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Consequently, since we have q > 3 and 2/p + 3/q < 1 by assumption we can choose m such that 
the constraint (j75p are matched. This yields that there exists 7 > 1 such that 



If Uq is sufficiently small, this yields that liuik^^oo Uk = 0. Since, we have from (|7ip . the Tcheby- 
chev inequality and the energy inequality (j68p that 

we can indeed make Uq arbitrarily small by taking M sufficiently large and thus limfc_^_|_oo Uk = 0. 
From Fatou's Lemma, we obtain that for every t G [0, 1] 

if{t,x) - M)+dx < 

and therefore that almost everywhere 

f{t,x)<M. 

By changing / into —f,we obtain in a similar way that / > —M almost everywhere and thus 
()67p is proven. To obtain the part of estimate (j65p involving i*" if T > 1 we can use a change 
of scale argument. Let us set Kf{T,X) = f{TT,VTX) fox K > to be chosen. Then we have 
KWfWifL^ = and f{T,X) solves the equation 

drf + u-Vx~f-^xf = yx-F 
where u is still divergence free and 

F{t,X) = ^F{tT,VtX). 



In particular, with the choice 
we get that Hi^HiPL.? = 1 and thus that 



1 3 \ 



L-L- = KWfU^L^ <MK = mVt^ VIIl^l.- 



This gives the part of the estimate (j65p involving F. 

Let us turn to the study of the second problem in (j66p in order to get the part of the estimate 
(j65p involving G. The estimate can be deduced from the previous one when qi < +00 which is the 
interesting case (when gi = 00, the estimate is a direct consequence of the Maximum principle). 
Indeed, if G G L^^L^i, we can write G = V -F with F £ L^W^'"?! C L'^L'^i by Sobolev embedding 
(where q^ = .) and moreover, we have 



Consequently, by using the estimate that we have already proven, we get that 



if 2/pi + 3/ql < 1. This gives the claimed estimate. 
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It remains to study the third problem in (j66p that is the problem with no source term but a 
nontrivial initial data. Again, we shall first prove that there exists M > such that for every 
/o G U with ||/o||l2 < 1, we have the estimate 

(76) sup||/(t)||Lco <M. 

i>l 

the standard energy estimate gives that 

(77) ll/llW + ||/||i.,.2<||/o||L- 

To improve this estimate, we shall also use the De Giorgi-Nash iteration method. We take a 
sequence M^, as previously, and we also choose a sequence of times = 1 — which tends to 1. 
The energy estimate for (/ — Mfc)_|_ yields that for every t, s with t >Tk > s, we have 

/ + 00 
||V(/ - M,)+)(r)||i, dT < IK/ - Mu)+{s)\\l, 
_ ,„ 

and hence, by integrating in s for Tfc_i < s < T^, we obtain that 

/■ + 00 

(78) Uk <{k + if / IK/ - Mk)+{s)\\l, ds 
with 

/■ + 00 

|2 



Uk = sup IK/ - Mk)+{t)\\i, + / ||V(/ - Mk)+)iT)\\i, dr. 



The aim is again to estimate the right-hand side of (j78l) by a power of Uk-i strictly greater than 
1. By using the same notations as previously, we get from ([72]) that 



(79) Uk<{k + If { ^^i±^) r t U- M,_,)J dtdx <ik + If ( ^±4^) ' Ul,. 



Since we have from ([77|) that Uq < \\{fo — {M — l))_^_\\i2 ^ ||/o|li2, Ui can be made arbitrarily 
small by taking M sufficiently large and hence we get from ([79j) that limfc__^+oo Uk = 0. This proves 
that f{t, x) < M for t > 1 and then we get (I76p by changing / into — /. We have thus proven that 
for every t > 1 the linear operator / i— )• /(t, •) is bounded from into L°° with norm smaller than 
M. Since by the standard maximum principle, it is also bounded from L°° to L°° with norm 1, we 
get by interpolation that it also maps L*" to L°° for every r > 2. To get the claimed estimate in 
()65p for t < 1, it suffices to use again a scaling argument. 
This ends the proof. 

□ 



Appendix B. Proof of Lemma [2771 
Proof. Set (p = J~^^h, then we have by definition and from Taylor formula 

[h{B)j]g = [ (^{x-y)g{y)[f{y)- f{x))dy 

= 11 g{y)^{x-y)-Vf{x + t{y-x))dydt, 

Jo JR-i 
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with = x(j){x). Let a,/3 G]0, 1[ with a + (3 = 1. Using Holder inequality and a change of 

variables we get with $t = 



[MD),/]5(x)| < C I {\9{y)mx-yT)[\Vf{x + t{y-x))mx-y)f)dydt 
Jo Jr'^ 

mx-y)\\g{y)\^)dyY [ [ I'^tiyW f{x - y)\hy 



< (ici>M5i^)"(x) r(i<i>ii*iv/p)''(x)dt. 

Jo 



Let pi,j)2 £ [1,00] such 



(80) api,/3p2 > 1 and - = — + —. 

P Pi P2 

Then by Holder inequality 

||[^(D),/]5||^, < \\m*\g\^\\l.r>, [' \\\^t\^\Vf\^\%,^dt. 

Jo 

We choose pi,p2,a and (3 such that 

111,111 

(81 am > 1, > 1, 1 + = - + and 1 + — = - + — 

api r am pp2 r pp 

then the classical convolution laws give 

\\[h{B)j]g\\^^ < mtrMi^w^ f\\LP [' \mldt 

Jo 

^0 

This last integral is finite despite that 

1 r 

(82) /3<-- 



3r - 1 



Now let us check that the set given by conditions (jSOj) . (|8T]) and (|82|) is not empty. First for the 
case r = 1 we choose pi = 'm,p2 = p,a = ^ and j3 = 1 — a. Let us now discuss the case r > 1. 
From (HI]) 

r 1 1 
a = -( ). 

r — 1 m pi' 

To get a G 1[ we must choose pi such that 

1111 

(83) < — < . 

p p pi rm 

The condition (3p> 1 is equivalent, by the use ofl + | = ^ + i + i, to 

,84) i-i<l. 

P rp pi 

The condition api > 1 is automatically satisfied from (|83|) since 

api > ram > r > 1. 
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The condition Pp2 > 1 is also a consequence of (j83j) and (|84j) . Indeed, from the value of a and 
I = ^ + this condition is equivalent to 

1 r 1 ^ 1 

p 2r — 1 p ~ pi 

This condition is weaker than (j84p . We can easily check that (j83p and (j84p are equivalent to 

1 1 ^ 1 ^ 1 

p rp ~ pi ~ rm 

The set of pi described by the above condition is nonempty if 

1 1 1 
+ — > 0. 

rm p rp 



Using the identity l + p = m~'"p~'"f' ^'^^^ satisfied under the condition p > r. The condition 

1111 
_ < + . 

pi 3 p p 



_i_ 1 _ 1 I 1 I 1 

(j82]) is equivalent to 

Now there is a compatibility between this condition and (j84p if 

3(1 - -) < p. 

r 

This ends the proof of Lemma 12.71 □ 
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